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Abstract 

To simplify a number of considerations in the weak field approximation, in- 
cluding the determination of the response of interferometric gravitational wave 
detectors, the "transverse traceless" (TT) gauge is often used. While the iden- 
tification of the corresponding gauge invariant part of the metric perturbations 
in the pure vacuum case is obvious, it is not widely known that the determina- 
tion and the use of the TT part is much more complicated and, in turn, much 
less rewarding when sources are involved. It is shown here that likewise the 
transverse part of the electric current in the Coulomb gauge within Maxwell's 
theory the sources of the TT gauge part of the metric perturbations become 
non-local. This, in practice, invokes the necessity of the use of more adequate 
projection operators then the ones applied, e.g, in the weak field limit, and 
in many post-Newtonian considerations. It is also pointed out that, whenever 
nonlinear effects are taken into account, some of the conclusions concerning 
the response of interferometric gravitational wave detectors may be influenced. 
In particular, attention is called on the possibility that gravitational radiation 
may produce an isotropic change of the spatial geometry. 
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1 Introduction 



General relativity is a metric theory of gravity which can also be formulated as being 
a gauge theory. This paper — besides pointing out various analogies between the use 
of the Coulomb gauge in Maxwell's theory and the TT gauge in the linearised Ein- 
stein's theory — is to point out the necessity of careful reinvestigation of the standard 
arguments applied in determining the response of interferometric gravitational wave 
(GW) detectors. The main motivation for the present work is rooted in the fact that 
sensitivity of the ground based GW detectors such as LIGO and Virgo has been im- 
proved significantly [UE1E1[Z]> and it is widely held that, if not earlier, then once the 
advanced detectors will be operating, detection of gravitational waves will become 
an everyday routine. Therefore, it is getting of obvious interest to have the best 
possible estimates for the astrophysical parameters of the associated GW sources. In 
this respect it is worth to be recalled that the optical observations of the change of 
the parameters of the orbital motion of potential GW sources agree — see, e.g., the re- 
ports on the Hulse- Taylor pulsar and similar type of binary systems [U El E] (see also 
section 6.2.3 of [12]) — , up to a very high precision, with the predictions of Einstein's 
theory. More precisely, the energy loss, which is signified by the observed change of 
the orbital parameters and which is assumed to be yielded by gravitational radiation, 
is in accordance with the predictions of Einstein's theory. To ensure the same order of 
precision in determining the astrophysical parameters based on the independent GW 
observations it is of crucial importance to be sure that not only the generation but the 
propagation of GW signals from the sources to our detectors is properly determined 
in the applied models. 

In this paper we intend to provide a simple enough discussion indicating some 
of the potential sources of imprecision in the current determination of the detector 
response for the arrival of GWs. In singling out an appropriate framework it turned 
out to be really useful to have the comprehensive paper by Flanagan and Hughes [8] at 
hand. Most of the arguments below are going to refer to the results formulated in the 
first part of this work. There is, however, one point to be mentioned here. As opposed 
to the ambitious plans manifested by the first part of [5] for certain reasons Flanagan 
and Hughes did not follow the path chosen there. For instance, after providing a very 
useful critical summary of the former conventional discussions in section 2.3 they 
returned to the orthodox arguments in spite of the fact that these are apparently 
inconsistent with the conclusion of the first part of their paper. In this respect, the 
main points of the arguments presented in this paper may be considered as natural 
continuation and completion of the work initiated by Flanagan and Hughes. 

As it was indicated above the use of the Coulomb gauge in Maxwell's theory shares 
several essential properties with that of the TT gauge in the weak field approximation — 
in particular, since the argument justifying the non- locality of the pertinent sources 
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in these gauges are completely parallel — in the rest of this section some of the most 
important related facts of Maxwell's theory are recalled. The gauge invariant Maxwell 
tensor F a p is given in terms of a vector potential A a as F a p = d a Ap — dpA ai and two 
vector potentials A a and A' a are known to be physically equivalent, i.e., they yield 
the same Maxwell tensor, if there exists a real function \ such that 

A' a = A a + d aX . (1) 

The field equations, whenever the gauge dependent vector potential A a satisfies the 
Lorentz gauge condition d a A a = 0, read as (see, e.g., [9| [TO]) 

□ A a = -4vr J a (2) 

where □ = —df + V 2 and J a stands for the electric four current vector. It is known 
that by choosing the real function \ appropriately the Lorentz gauge condition can 
always be guaranteed to hold. Recall also that we still have the freedom of applying 
a restricted gauge transformation of the form ([1]) provided that the generator \ is 
subject to the equation □ \ = since then the Lorentz gauge condition remains 
intact. 

Is is also well-known that gauge independent quantities can be built up from the 
vector potential. The pertinent gauge is frequently referred as "Coulomb", "radia- 
tion" or "transverse" gauge and it can be introduced as follows. Start by picking 
out an inertial reference system, (t, x), of the underlying Minkowski spacetime. Then 
a vector potential A a may be decomposed^ as A a = (— 0, Aj). The spatial part A^ 
of A a can be split up into 'transversal' and 'longitudinal' part as A^ = Aj + dtp, 
where Aj is such that d l Aj = 0. This decomposition is unique if in addition the 
potential p is guaranteed to tend to zero while r — > oo as then the elliptic equation 
V 2 p = d l Ai possesses a unique solution. It is straightforward to see that once a gauge 
transformation (pQ) is applied the variables <fi and </? will be changed. It is well-known, 
however, that their combination $ = + d t f, along with Aj, is gauge invariant [9]. 
The field equation ([2]), pertinent for $ and Aj, reads then as 

V 2 $ = -47rp (3) 

nA? = -4*j? (4) 

where the decomposition J a = (— p, Jj) has been used. Notice that the source for the 
transverse part of the vector potential Aj, 

j? = Ji--}-m*) (5) 



lr The minus sign in front of the scalar potential comes from the fact the for a vector potential 
A a . with Aq = —(f>, entered into the formalism much later then the Coulomb potential. 
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extends over all space even if the spatial part Jj of J a is localised [9]. This is a 
direct consequence of the fact that $ is subject to the Poisson equation (jSJ), i.e., $ is 
non-local. Notice that $ is time dependent even though it does not time evolve as a 
wave. 

An immediate consequence of the non-locality of the transverse part of spatial 
vector fields is that, once sources are involved, one has to be careful in determining, 
e.g., the transverse part Aj of the spatial part of vector potential Ai by making 
use of a projection operator. A projection operator!! Pi of this type is formally 
defined referring to the inverse Laplace operator, ^2, as Pi = Si — di^&>. However, 
the precise form of this projection operator — see, e.g., the discussion on page 242 
of Jackson's book [9J — , making the non-locality of the involved fields completely 
transparent, assigns to a spatial vector Vi its transverse part = PiVj, as 

W = w + ±pf^l**. (6) 

where the relation V x | x _ x /| = ^( x ~~ x has implicitly been used. 

It is straightforward to verify that whenever the spatial vector Vi possesses the 
form of a plane wave solution, i.e., V, = V® ■ cos(kx — ut + ip ), with constant 
amplitude V z ° and phase ipo, and with uo = |k|, then Pi can be given as 

Pi j = <V - rnn j , (7) 

where the spatial unit vector rij is given as rti = ki/u. It is important to emphasise 
that the projection operator Pi , given by ([6]) for the generic case, reduces to the form 
of ([7]) if and only if Vi is given as a linear superposition of plane wave solutions such 
that all the spatial wavenumber vectors are parallel. Accordingly, the application of 
the projection operator Pi (j7|) does not yield the TT-part of Vi besides this excep- 
tional case. It is worth to be mentioned that it does not even do the job for slightly 
more general solutions to the sourceless wave equation, □ Vi = 0. All these obser- 
vations imply then that whenever sources are involved the only adequate projection 
operator must possess the form of (jSJ). 

2 The weak field approximation of GR 

The weak field approximation of general relativity is believed to be adequate in de- 
scribing weak gravitational effects. In such a case the metric g a p of the spacetime 

2 All the index raising and lowering are meant to be done by either of the fixed background 
metrics r\ a f$ or Sij of the Minkowski spacetime or the Euclidean space, respectively. Moreover, the 
Einstein's summation convention is used only for identical upper and lower indices. 
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is supposed to be close to the flat metric r\ a p of the Minkowski spacetime. More 
precisely, it is assumed that Minkowski-type coordinate systems exist such that 

g a p = Va(3 + h a p (8) 

and that \h a p\ <C 1 • 

As a direct consequence of the generic diffeomorphism invariance of Einstein's 
theory two linear perturbations h a p and h' a/3 of the flat Minkowski spacetime are 
considered equivalent, whenever they are related as 

Kp = h *P + d a ^p + dp£ a , (9) 

where £ Q denotes some infinitesimal vector field determining the coordinate transfor- 
mation 

x a^ x >* = x *-£\ (10) 

Notice that in (JHj) h a p and £ a play the same role as the vector potential A a and the 
function x do in the Maxwell case. 

The linearised Einstein equations can then be shown to take — in terms of the trace 
reversed, 

Kp = h aP - - r] a(3 /i 7 7 , (11) 

of h aj 3 — the simple form 

UKp = -lQnT aP (12) 

provided that h a p satisfies the Lorentz gauge condition 

d a h a p = 0. (13) 

It is well-known but worth to be mentioned that there always exist coordinate trans- 
formations of the form x' a = x a — £ a such that ([13"]) holds in the new gauge. Moreover, 
the pertinent gauge is not unique since further restricted coordinate transformations 
with £ a subject to the wave equation 

DC = (14) 



may still be applied as they leave the Lorentz gauge condition f|T3|) intact. 

The solution to the inhomogeneous equation (I12p . given in terms of the retarded 
Green function, read as 



T a p(t - 




x — 


x' 


,x') 




x 


-x' 





h a p(t, x) = 4 / ^ PV ", l - ""-' rfV . (15) 



In virtue of (TT2I) all the components of h a p possess radiative degrees of freedom which 
by many authors (see, e.g., [SI ED]) is considered to be an "unfortunate consequence" 
of the applied gauge. A more adequate objection could be that the components of 
h a p are not gauge invariant thus they cannot be directly applied in determining the 
response of our GW detectors. 
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3 The "radiation" or TT gauge 

In virtue of the criticism recalled above, more importantly, because of the obvious 
need for a correct derivation of the response of our ground based laser interferometric 
detectors like LIGO and Virgo to the arrival of a GW signal, it is important to know 
whether the true radiative physical degrees of freedom can always be separated in the 
weak field approximation. 

It has been known for long that the gauge independent expressions can be built 
up from the components of h a g. In the following short review of the pertinent results 
we shall follow the discussion of [8] unless otherwise indicated. 

Consider first a "1 + 3" decomposition 



h 



a/3 





hu 


hu 


h^ 



(16) 



of h a/ 3 based on the use of a Minkowski type coordinate system, (t, x), where time- 
time, time-space and space-space parts are given in terms of the variables 0, 7, £j, 
hJ, T and A as 

hu = 20 (17) 
thi = Pi + do (18) 

h n = ^ T + \ H $ij + + \9idj - l^V 2 ^j A , (19) 

where H = 5 l ^hij denotes the three-dimensional trace which is related to h = h a a as 
h = H — 20. The variables 7, £j, A — and, in turn, $ and hj? as well — can be seenEl 
to be uniquely determined once the relations 

d^i = , d% = , d^Jj = , (20) 

along with the boundary, or fall off, conditions 

7^0, £i -> 0, A -> 0, V 2 A -> while r -> 00 , (21) 

are imposed. Note that, in virtue of fTl9l and (120]) the contraction S^hJ? vanishes, 
which along with the last relation of (120]) . implies that hj? is TT. 

As the components of h a p themselves are not gauge invariant the variables 0, 7, A, 
H, (3i and £« are not gauge invariant either. However, the combinations 

$ee-0 + %-^ 2 A (22) 
9 = i (H - V 2 A) (23) 
Hi = A - , (24) 



3 For more details the reader may look up the pertinent part of the argument applied for the 
analogous decomposition of the energy-momentum tensor below. 
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along with the 3x3 matrix tiff, can be shown to be gauge invariant. 



3.1 The decomposition of the energy- momentum tensor 

In order to be able to determine the evolution equations for the above introduced 
gauge invariant expressions we shall need an analogous decomposition of the energy- 
momentum tensor. 

Before providing this decomposition recall first that whenever matter fields are 
involved the Einstein's and matter field equations have to be solved simultaneously. 
Now, partly to simplify our argument, and also to avoid the associated considerable 
technical difficulties, without choosing any concrete field equations we shall assume 
that the field values, along with the components of the energy- momentum tensor, T a p, 
are determined by some unspecified field equations, governing the time evolution of 
the sources of the gravitational waves. Once we have the energy-momentum tensor, 
T a/ 3, a decomposition, completely analogous to the one applied above for h a p, can be 
provided as follows. 

Start by a "1 + 3" splitting of T a p 



T, 



a/3 





T tl ' 


T u 


Tij 



(25) 



and by defining the variables p, Si, S, a^, Oi and a via the relations 



Tu = P (26) 
T ti = S t + (27) 

T^ = <Tij + PSij + d (i a j} + {did, - ^yVj o , (28) 

where P = ^d^T^ . As above S, a and Oi get to be uniquely determined once the 
relations 

&Si = , &(n = , &<Tij = , (29) 
along with the boundary, or fall off, conditions 

S -> 0, Oi -> 0, a -> 0, VV -> while r -> oo , (30) 

are imposed. Note that as above, in virtue of (|28[) and (I2U1) the contraction S^a^ 
vanishes, which along with the last relation of (129]) . implies that &ij is TT. 

The uniqueness of the above decomposition can be seen as follows. First, the 
^-divergence of (1271) yields S7 2 S = d l T ti , which has a unique solution by the above 
boundary condition. Once S is known Si gets to be uniquely determined by the rela- 
tion Si = Tu — diS. Concerning the uniqueness of a, take now the d l &>- "divergence" 
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of which yields the Poisson equation V 2 V 2 cr = § [d l d j Tij — V 2 P], and which 
has a unique solution for V 2 cr, and, in turn, in virtue of f[3"Ul) . cr becomes uniquely 
determined, as well. Once V 2 a is known <jj gets also to be uniquely determined 
by the Poisson equation V 2 ^ = 2 [&>Tij — diP] — | <9i V 2 cr, which is yielded by the 
<9 J -divergence of f[2"8"]) . along with the fall off condition Oi — > while r — > oo. 

As the energy momentum tensor is supposed to be known it is straightforward to 
see that having a and Oi determined, cr^- gets also fixed as 

<Tij = Tij - P5 i:j - d(j<Tfi - (didj - ^yV 2 J a . (31) 

In [8] the authors claim that <7y, along with some of the other variables, can be chosen 
freely It should be noted that this can be done only if the energy-momentum tensor 
is not specified. 

There is an even more important additional point to be mentioned here. As the 
variables a and <jj satisfy Poisson type equations they are non-local. In consequence 
of this fact and the above relation (j31j) the TT-part, cr^ , of cannot be local either 
regardless whether the energy-momentum tensor, T a p, of the matter sources is of 
compact support or not. 

Since this non-locality is in certain extent inconvenient it could be tempting to 
argue that although the variables S, a and cr, were shown to be subject to some very 
complicated elliptic equations they might be completely negligible. In this respect it 
is useful to take into account the conservation law d^T^ = which reads as [8] 

V 2 S = d t p (32) 
V 2 a = - 3 -P+ 3 -d t S (33) 
V 2 a i = 2d t S i . (34) 

These relations — which had a completely different role in the discussion of [S]- 
provide the following alternative characterisation of the variables S, a and <7j. They 
make it immediately transparent that the time derivative of the energy density and 
that of the impulse — both of these quantities are supposed to provide significant 
contribution to the gravitational wave production — are the sources for S and a"j, 
respectively. Thereby, the more intensive is the considered GW source the more 
significant the quantities S, a and cr, become. 

4 The linearised Einstein's equations 

Now we are prepared to present the explicite form of Einstein's equations relevant 
for the above introduced gauge invariant quantities in the investigated weak field 
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approximation. These equations can be given as [8] 



V 2 6 = -8ttp (35) 

V 2 $ = An(p + 3P-3d t S) (36) 

V 2 H, = -lQnS t (37) 

□ \5 r = -167r«7y. (38) 



The above equations justify the conventional assertion that only the TT part of 
the metric perturbation satisfies wave equation while all the other gauge invariant 
expressions, although they are time dependent, do not time evolve as waves since 
they are subject to Poisson equations. Accordingly, it is frequently said that only 
the "non-radiative" physical degrees of freedom are tied to the matter sources. What 
is even more surprising is that conclusions of the following type are drawn based 
on the above set of equations: Since the sources are at enormous distance from the 
Earth, in virtue of GW signals can basically be considered as being sourceless 

and possessing the same type of properties as if they were GW signals in the pure 
vacuum case. As opposed to this, we would like to emphasise that according to the 
conclusion of the previous section the source term cry in fl38l) is non-zero, and it is 
non-local either even though T a p is completely localised. 

As an immediate consequence of this non-locality let us make a comment regarding 
the "conventional" way of determining the TT part hj^ of a solution h a p to the 
evolution equation (jl2j) . It is usually assumed in the weak field approximation (see, 
e.g, Section 4.1 of [8] or section 3.1 of [12]) and, for some surprise, also in the post- 
Newtonian framework (see, e.g., the sentence involving Eq. (2.2) of |13| . section 5.3.4 
of [12] or the first paragraph on page 19 in |14j). that the TT part of may be 
determined in terms of the projection tensor, 

Kj kl = P k P 3 l - l -P l3 P kl , (39) 

as hJ 3 T = Aij kl hki, where the "elementary projection operator" Pi is supposed to 
possess the form Pi = Si — n^nK We would like to emphasise here that, as it follows 
from the argument outlined at the end of section[TJ this form of Pi cannot adequately 
be applied even to the superposition of plane wave solutions to the sourceless wave 
equation unless all the spatial wave number vectors are parallel. Since the evolution 
equation ([12]) for h a p must have sources in astrophysical situations, as actually we wish 
to observe these sources, the projection operator Ki k , when it is expressed in term of 
the "elementary projection operator" via (1391) — as opposed to the generic resolution 
applied in various calculations (see, e.g., [HI d21 EH EH]) — has to be constructed by 
making use of ([6]) instead of applying ([7]). The corresponding projection operator — 
the complexity of which is expected to reflect all the technical difficulties related to 
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the non-locality of <7y — will produce the adequate TT-part of h\j. In virtue of these 
observations there is an obvious need for the reinvestigation of the procedures yielding 
the wave forms in the asymptotic region by applying the correct projection operator 
Ajj M , e.g, in the post-Newtonian framework. 



5 Further implications of non-localities 

In proceeding, let us recall now, that in many of the arguments, aiming to determine 
the response of laser interferometric detectors to the arrival of a GW signal, the 
calculations end up (see, e.g., Eq. (3.10) of [8], or Eq. (1.93) of [12]) with the variant 
of the geodesic deviation equation 

where L l {t) = L % + SL l (t), with 5L l {t) -C L % and i,j = 1,2, is supposed to denote 
the coordinates of mirrors at the end of the arms in the proper detector frame. 

Then, in the linearised theory, assuming that no sources are present, the relation 
Rujt = — \^t^i7 1 a l° n g with the assumption that both 8U(t) and d t (5L l (t)) vanish 
at t — 0, is applied to derive the familiar "gauge independent" relation 

6Li{t) = \l^li. (41) 

However, as it has already been emphasised above, since we do want to make 
astrophysical observations, the presence of the sources has to be taken into account. In 
this more realistic situation, as opposed to the pure vacuum case, the gauge invariant 
"tidal force components" of the Riemann tensor read as (see, e.g., [8]) 

Rm = + d ^ + d A^i) ~ \ %® Stj ■ (42) 

It is well-known that in the linearised theory, if one takes into account the conservation 
of the stress energy tensor, d a T a p = 0, the last three terms can be shown to fall off 
like 1/r 2 or faster [8J0 Thus, the only term with 1/r fall off is the first term on 
the right hand side of (142 1) . Thereby, within the linearised theory regardless whether 
sources are present relation (I4ip determine the response of our detectors. 

Up to now only the linearised theory has been considered. However, it has been 
known for long (see, e.g., the discussion in section 4.4. of |10j ) that it has serious 



4 In virtue of (|35|) - (|37p <&, 8 and Ej fall of like 1/r. However, the coefficients of the 1/r parts of 
these quantities can be given in terms the conserved mass M = J pdx and the conserved linear 
momentum Pi = J Si dx 3 [5] . 
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limitations, since, even to have a consistent Newtonian limit one must go beyond the 
linear approximation. This is justified by the fact that in the linear theory, in virtue 
of the conservation law cTT^p = 0, e.g., the members of a binary system — instead of 
orbiting around each other — have to follow timelike geodesies, i.e., straight lines, of 
the Minkowski spacetime. As opposed to this a nearly Newtonian type of orbiting is 
produced by the binary if the terms higher order in h a p are restored in the Einstein 
tensor G a p. Therefore, in a physically consistent description of gravitational wave 
generation processes the back-reaction has to be taken into account. This can be done, 
while preserving the simplicity of the basic equations of the linear approximation, by 
putting all the non-linear terms of the Einstein tensor to the energy-momentum tensor 
side, or, more precisely, by replacing in ffT2]) the energy-momentum tensor T Q/ g by the 
sum T a/ 3 + t a p, where t a p = — ^ ln °G Q/ g, and {nl) G a p consists of all the higher order 
terms in the Einstein-tensor. Note that the assumption, \h a p\ <C 1, requiring the 
perturbations to be sufficiently small, could also be relaxed then. In particular, once 
T a/3 is replaced by the sum T a p + t a p the yielded equations become suitable to describe 
the evolution of intrinsically strong GW sources which cannot be done properly in the 
linearised theory. Note also that the conservation law c^T^ = 0, which is responsible 
for the above mentioned defects, gets to be replaced by the more adequate relation 
d a (T a/ 3 + t a p) = 0. Once the replacement T a p — > T a p + t a p, along with a senseful 
redefinition of the quantities p, Si, S, P, a^, <r i? has been done then all the previously 
derived equations, (j52"|) - (j58|) . can be seen to preserve their forms. 

Now, by making use of the above introduced nonlinear setup, we intend to provide 
a plausibility argument suggesting that the response of our detectors is going to be 
affected by back-reaction. Before presenting our argument we would like to emphasise 
that, whenever nonlinearities are taken into account but sufficiently far from the 
sources \h a p\ <C 1 holds, it seems to be reasonable to assume that the tidal forces 
can still be given by (|42p . with the distinction that now and hj^ refer to 

the redefined quantities. We shall use this assumption below. Note also that in 
consequence of the hidden nonlinearities no attempt is made to go beyond providing 
a plausibility argument, i.e., no quantitative estimates are derived. 

In proceeding let us revisit the fall off properties of the last three terms on the right 
hand side of (I42p . Note first that since the middle two terms in (1421) contain spatial 
derivatives, it is straightforward to verify that they decay faster than the terms on 
the sides. In particular, since both $ and Hj fall off like 1/r the relations didj§ ~ r -3 
and dtduB,j\ ~ r~ 2 can be seen to hold. Therefore the main issue is whether, besides 
hj? ', the last term on the right hand side of (]42[) may also have a 1/r fall off in the 
nonlinear case. In this respect the following simple example provides some important 
clues. 

Assume that in our spacetime we have nothing else but a localised GW source 
which produces a single short lasting GW burst. Consider now an observer that is 
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asymptotically far from the source in the distance. For simplicity let us represent 
the world-sheet of the GW burst, as it is travelling from the source towards infinity, 
by a null shell. The observation occurs where the world-line of the observer meets 
this null shell. In advance to the observation the 1/r part of G (as measured by the 
observer) is proportional to the ADM mass. However, if back-reaction is taken into 
account, regardless how tiny is the energy carried by the GW on the future side of 
the null shell, the coefficient of the 1/r part of G will be smaller. This is so because 
the mass felt by the observer will be smaller than the ADM mass as the GW, even 
after its detection, goes on, carrying the energy released by the source, towards null 
infinity. Accordingly, the coefficient of the 1/r part of G will vary in time in spite of 
the fact that the total ADM mass is conserved. 

The above example suggests that nonlinearities manifest themselves in the follow- 
ing simple way. Whenever T a/3 is replaced with T a p + t a p — in consequence of the fact 
that tap extends beyond the observer and it is time dependent — the coefficient of the 
1/r part of Q, at the location of the observer, vary with time. If the observer is not 
too far from the sources then this time variation might be oscillatory although it is 
expected to be monotonic, as in the above example, if the observer is asymptotically 
far from the sources. 

Thereby, in the nonlinear regime it seems to be reasonable to assume that second 
time derivative of both G and hj? fall off like 1/r. Accordingly, in virtue of ()42|) . for 
the variation of the coordinates of the end mirrors the relation 

8Li(t) « i [hj^ + QSij] L{ (43) 

applies. What is even more remarkable is that the effect of the second term on the 
right hand side of (|43i) on the arms of the laser interferometric detectors is nothing 
but a possibly tiny but isotropic change of the arm lengths. 

Of course, without carrying out further quantitative investigations there is no 
way to argue that this effect is important. Thus, it is of obvious interest to know 
whether this effect is large enough, and whether it could be detected by the current 
arrangements of our ground based laser interferometric detectors. Unfortunately, 
the answer to the latter question is no since the LIGO- Virgo type detectors in their 
present form are sensitive only to the relative variation of the arm lengths. We hope 
that by a suitable modification of the applied detection schemas, e.g., by adopting 
some of the ideas proposed in [TS], it may be possible to measure the variation of 
the arm lengths separately. If this can be done an isotropic change of the spatial 
geometry, simultaneous to the arrival of a gravitational wave train, could hopefully 
be detected. 
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6 Final Remarks 



In this paper some of the peculiarities of the TT gauge in the weak field approximation 
were investigated. The results found indicate that there have been several assump- 
tions applied in determining the TT part of the metric perturbations which are in 
the air. In particular, the results found have the following non-trivial consequences. 

First, it is pointed out that, whenever sources are involved, in determining the 
TT part of the metric perturbations a new approach — significantly different from the 
currently applied one — , which is taking into account the non-locality of the sources 
pertinent for the gauge independent TT variables is needed. It is also indicated that 
the associated improvements will affect the asymptotic wave forms not only in the 
linearised theory but in the post-Newtonian framework, as well. 

Second, the determination of the response of interferometric gravitational wave 
detectors may be influenced considerably if back-reaction is taken into account. A 
plausibility argument was provided justifying that the arrival of a GW train yields 
an additional isotropic change of the arm lengths. We would like to emphasise that 
there is an immediate consequence of this effect which may affect the estimates con- 
cerning the detectability of gravitational wave signals by our current detectors. This 
is related to the possibility that the energy released by the astrophysical sources may 
not completely be transferred into the pure radiative degrees of freedom but some 
part of it could be used to produce the isotropic change of the spatial geometry in 
the distance. This variation of the volume may decrease the current estimates of GW 
amplitudes given in terms of hj^. It is also indicated that the anticipated isotropic 
change may be observed by suitably modified versions of the currently applied laser 
interferometric detectors. 

It is worth to be mentioned that in most of the investigations of GW productions, 
whenever there is an attempt to take into account the nonlinear back-reactions, cu- 
riously enough, the nonlinearities themselves are left out almost immediately from 
the discussions. For instance, in determining the quadrupole tensor, which, in the 
nonlinear case, should read as 



it is usually assumed (see, e.g., the bottom of p. 87 in [10], or section 4.2 in [8]) that 
the contribution of t 00 can be neglected as the relation t 00 -C T 00 holds at the location 
of the source. However, in accordance with the last remark of Wald on the top of p. 88 
of [TU], we would like to emphasise that a higher level of clarity and rigour should 
be involved here. Note that the main results of the present paper support these 
necessities simply because t 00 is known to be global, and, far from the sources, where 
t 00 may be small but positive, its contribution could also be significantly amplified 




(44) 
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by the factor x a x 13 in (jUJ). Therefore, it may happen that whenever t 0Q is sufficiently 
anisotropic its contribution to the quadrupole moment tensor will become significant. 
This, in turn, may yield a considerable deformation of the emitted wave forms while 
travelling from the sources to the observers. If this turns out to be the case we cannot 
avoid a careful revision of the currently applied template banks. 

It is also of obvious interest to know what might be the relation between our 
findings and "Christodoulou's nonlinear memory effect" [T5]. In this respect we would 
like to mention that, because of the significant differences of the mathematical setup 
applied in [16] and in this paper, it is not obvious at all to derive a meaningful 
relation especially because in [16] no explicite expression is given that could represent 
the displacement of the mirrors. In particular, as a consequence of these differences 
while in this paper the nonlinear effects were shown to be isotropic no such conclusion 
was derived in [IB]. Nevertheless, it is worth to be mentioned here that according 
to the rough estimate provided by Christodoulou the nonlinear effects may be of the 
same order as the linear ones. 

It might be tempting to consider further consequences of the indicated isotropic 
change of the spatial geometry. Indeed, the pertinent implications might be far- 
reaching if at certain parts of the universe the monotonous increasing of the function 
G could be verified. Nevertheless, we would like to emphasise again that in this respect 
the present paper is far from being conclusive as neither reliable estimates concerning 
the fraction of the energy converted into the change of has been derived nor its 
monotonicity has been studied. To do so further analytic and numerical investigations 
are needed. 
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